Two-stream Instability in Pulsar Magnetospheres 



O 
O 

^ ■ 

Oh; 



> ■ 
<N : 
o 

o 
o 

Oh: 

6 



V.V. Usov 

Department of Condensed Matter Physics, Weizmann Institute of Science, Rehovot 76100, Israel 



Abstract. Creation of electron-positron pairs near the 
pulsar surface and the parameters of plasma in pulsar 
magnetospheres are discussed. It is argued that the pair 
creation process is nonstationary, and the pair plasma that 
flows out from the pulsar environment is strongly non- 
homogeneous and gathers into separate clouds. Plasma 
instabilities in the outflowing plasma are reviewed. The 
two-stream instability that develops due to strong nonho- 
mogeneity of the outflowing plasma is the most plausible 
reason for the generation of coherent radio emission of 
pulsars. The development of the two-stream instability in 
pulsar magnetospheres is considered. 



1. Introduction 

Current models for the generation of coherent radio emis- 
sion of pulsars require the development of plasma insta- 
bilities in the pulsar magnetospheres, the most studied 
of which is the two-stream instability (for a review, see 
Arons 1981a,b; Melrose 1981, 1993, 1995; Usov 1981; As- 
seo 1996). The two-stream instability in the pulsar magne- 
tospheres has been considered in many papers (e.g., Rud- 
erman & Sutherland 1975; Benford & Buschauer 1977 
Buschauer & Benford 1977; Cheng & Ruderman 1977 
Arons 1981; Usov 1987; Ursov & Usov 1988; Asseo 1993 
Asseo & Melikidze 1998; Melikidze, Gil & Pataraya 2000). 

The two-stream instability can arise when two plasma 
flows travel through each other: 



vi, m 



V2, "2 



If the density of one flow is much higher than the density 
of the other (for example, ni 3> n2), the low-density flow 
(component 2) calls the beam while the high-density flow 
(component 1) calls the plasma. In the magnetospheres 
of pulsars, both the plasma-beam interaction and the in- 
teraction of two plasmas with more or less equal densi- 
ties may occur (see below). The two-stream instability 
may lead to formation of plasma bunches that generate 
the radio emission via curvature mechanism (Ruderman 



& Sutherland 1975; Cheng & Ruderman 1980; Ochelkov & 
Usov 1984). Besides, the radio emission of pulsars may be 
generated directly in the process of the instability devel- 
opment (Asseo, Pellat & Rosado 1980; Asseo, Pellat & Sol 
1983). Before to discuss the two-stream instability in the 
pulsar plasma and the generation of coherent radio emis- 
sion of pulsars, we review both the physical processes that 
are responsible for the filling of pulsar magnetospheres 
with plasma and the plasma parameters. 



2. Physical processes and parameters of plasma in 
pulsar magnetospheres 

A common point of all acceptable models of pulsars is that 
the spin-down power of strongly magnetized neutron stars 
is the energy source of non-thermal emission of pulsars. 
The rotational energy of the neutron star is transformed 
into the pulsar emission by a long sequence of processes. 

2.1. Physical processes in pulsar magnetospheres 

Strong electric fields are generated in the magnetospheres 
of rotating magnetized neutron stars (for a review, see 
Michel 1991). The component of the electric field E|| = 
(E • B)B/|Bp along the magnetic field B is non-zero. Pri- 
mary particles are accelerated by this electric field to ul- 
trarelativistic energies and generate 7-rays. Some of these 
7-rays are absorbed by creating secondary e^e~ pairs. 
The created pairs screen the electric field E|| in the pul- 
sar magnetosphere everywhere except for compact regions. 
The compact regions where E|| is unscreened are called 
gaps. These gaps are "engines" that are responsible for 
the non-thermal radiation of pulsars. The gaps are located 
either near the magnetic poles of pulsars (polar gaps) or 
near their light cylinders (outer gaps). Outer gaps may 
act as a generator of non-thermal radiation of pulsars 
only if the period of the pulsar rotation is small enough, 
P < Per — a few X 0.1 s. For typical pulsars with P > Per, 
the polar gap model has no an alternative. In the polar 
gap model the sequence of processes that leads to the gen- 
eration of non-thermal emission of pulsars is the following: 
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Physical processes in pulsar magnetospheres 



generation of strong 
electric fields, 



Eii = 



(E-B) 

I B I 



#0 



generation of 
electromagnetic 
oscillations 



coherent radio 
emission 



acceleration 
of primary 
particles 



development 
of plasma 
instabilities 



generation 
of 
y- rays 



creation of secondary 

e"*", e" plasma, 
Y+ B— ►e"'" + e" + B 



non-thermal high-friquency 
(IR, optical. X-ray and y-ray) 
emission 



The density of primary particles accelerated in the 
polar gap is about the Goldreich-Julian density ncj = 
jn • B|/27rce, where ft is the angular velocity of the pul- 
sar rotation. The value of ncj is determined so that the 
electric field E^^ in the outflowing plasma is screened com- 
pletely if the charge density is equal to encj- Therefore, 
the accelerating field -Ey arises from deviations from this 
density. Such a deviation may be because of (1) the in- 
ertia of particles (Michel 1974), (2) the curvature of the 
magnetic field lines (Arons 1981), the General Relativity 
effects (Muslimov & Tsygan 1992, Beskin 1992), and the 
binding of particles within the solid surface of a strongly 
magnetized neutron star (Ruderman & Sutherland 1975; 
Cheng & Ruderman 1980; Usov & Melrose 1995, 1996). 
When both the neutron star surface is cold enough and the 
surface magnetic field is strong enough, there is no emis- 
sion of particles from the neutron star surface, and the 
binding of particles is responsible for the field in the 
polar gaps (e.g., Usov & Melrose 1996). In the case when 
particles flow freely from the stellar surface, the field E^^ 
is mainly due to the General Relativity effects (i.g., Hard- 
ing & Muslimov 1998). The second case is applicable to 
typical pulsars while the former case may be applied to 
pulsars with very strong magnetic fields at their surface. 



10^^ G or more. 



Ultrarelativistic primary particles in the process of 
their outflow move practically along the magnetic field 
lines and generate 7-rays via both curvature radiation 
and magnetic Compton scattering. These 7-rays are ab- 
sorbed in strong magnetic fields of pulsar magnetospheres 
and create e+e" pairs. Created pairs very quickly lose the 
momentum component transverse to the magnetic field 



due to synchrotron losses, and their distribution becomes 
one-dimensional, i.e., v || B. Hence, to develop the the- 
ory of coherent radio emission of pulsars it is necessary 
to investigate both the parameters of ultrarelatistic one- 
dimensional plasma outfiowing from the polar cap regions 
and instabilities of this plasma. 

2.2. Parameters of pulsar plasma 

Electrons and positrons that flow away from the neu- 
tron star vicinity with relativistic speeds may be di- 
vided into two components: a secondary e~^e~ plasma and 
an extremely high-energy beam of primary particles (see 
Fig. 1). Particles of the beam may be either electrons or 
positrons, depending on the direction of the electric field 
E|| in the polar gap. For the e+e" plasma (denoted by p) 
and the ultrarelativistic beam (denoted by b) the typical 
values of the Lorentz factor F and the density n are, re- 
spectively, (Ruderman & Sutherland 1975; Arons 1981a, b, 
1984; Machabefi & Usov 1989; Michel 1991) 



Fp ~ (2 - 3)Rc/R ~ a few X (1 
Tb eAip/meC^ ~ a few x 10^ , 



10^), 



Up ~ {Tb/2Tp)nb , 

Ub n^j = QB/{2Trce) ~ 10^\R/rf cm"^ , 



(1) 
(2) 
(3) 
(4) 



We have adopted here the following characteristic values 

for the pulsar parameters: ~ 10 s"-'^; B ~ 10"'^^(i?/r)^ G 
is the magnetic field in the magnetosphere; Rc is the cur- 
vature radius of the magnetic field lines in the polar cap; 
R ~ 10^ cm is the neutron star radius; and Atp is the 
potential drop across the polar gap. 
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Fig. 1. The energy spectrum of the basic components of the plasma particles ejected from the polar gap region into 
the pulsar magnetosphere. The spectra of secondary electrons (dashed line) and secondary positrons (solid line) are 
shifted relatively each other. 



2.3. Nonstationarity of the plasma outflow 

It was suggested long time ago that the process of pair 
creation near the pulsar surface is strongly nonstationary 
(Sturrock 1971; Alber et al. 1975; Ruderman and Suther- 
land 1975). Such a nonstationarity was discussed in detail 
in the outstanding paper by Ruderman and Sutherland 
(1975) where the first self-consistent model of pulsar radi- 
ation was presented. In the Ruderman-Sutherland model, 
the creation of pairs occurs in narrow spark discharges, 
separated from one another by a distance of the order of 
the polar gap height H (characteristic value of if is ~ 10^ 
cm). The secondary e+e~ plasma, produced by absorption 
of 7-rays in a spark discharge, gathers into an outflowing 
jet that has the discharge at its base. Such a discharge is 
strongly nonstationary on a time scale of ~ 30 (if/ c) (Rud- 
erman fc Sutherland 1975; Cheng fc Ruderman 1980). The 
discharge nonstationarity results in a strong modulation 
of the outflowing pair plasma along the jet. In fact, these 
jets consist of plasma clouds spaced by Li ~ 30H ~ 0.3R. 

However, it is not only the nonstationarity of dis- 
charges that can lead to inhomogcncity of the outflow- 
ing plasma along the fixed magnetic field lines. Another 
reason for such a inhomogeneity is the displacement of dis- 
charges caiiscd by the curvature of the magnetic field lines 
(Ruderman & Sutherland 1975). The displacement is di- 
rected toward the magnetic axis of the pulsar. The speed 
with which a spark discharge moves across the polar cap is 
estimated (Ruderman & Sutherland 1975; Cheng & Rud- 



erman 1980) as Vd — cH/2Rc. As already noted above, the 
average distance between two discharges is ^ _ff . During 
the time Td — H/vd — 2Rc/c, a discharge is displaced by 
the distance of ^ iJ and takes a position close to the ini- 
tial position that of a neighboring discharge, which during 
the same time is also displaced hy ^ H toward the pul- 
sar magnetic axis. As a result of these displacements of 
discharges, the plasma outfiowing along any fixed mag- 
netic field lines becomes modulated with a characteristic 
length L2 — cTd — 2Rc. Since, according to Ruderman and 
Sutherland (1975) and Beskin (1982), in the polar gap the 
discharges are separated by intervals within which the pro- 
cess of pair creation is suppressed, the plasma density in 
between the clouds falls down almost to zero. Therefore, 
the outflowing plasma gathers into separate clouds spaced 
by L2 — 2Rc along the direction of the plasma outflow. In 
turn, these large clouds consist of small clouds with the 
lengths and the spaces between them about Li ~ 0.3R. 

The idea that the process of pair creation near the 
polar caps of pulsars is nonstationary has got recently the 
strong support (e.g., Muslimov & Harding 1997; Harding 
& Muslimov 1998). In these papers, the model of the polar 
gaps of typical pulsars which is rather close to reality is 
developed. In this model, it was shown that both stable 
acceleration of primary particles and stable creation of 
pairs are not possible. The characteristic length of the 
outflowing plasma modulation is expected more or less 
comparable with Li. 
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3. Possible instabilities in pulsar plasma 

Development of instabilities in a magnetized plasma when 
transverse electromagnetic waves (i-waves) are generated 
depends on the ratio Wp/ws, where 



= / 5.64 X lO^nV^ , (5) 

y Tile 

eB 

ujb = ^ 1.76 X lO^B s"^ , (6) 

UleC 

are the Langmuir frequency of the plasma and the cy- 
clotron frequency, respectively. At the pulsar surface for 
the typical parameters of both the pulsar plasma and 
the strength of the magnetic field {up^ ~ 10^'* cm~^ 
and ~ 10^^ G), from equations (5) and (6) we have 

Since the difference between the phase velocity {vph) 
of waves and the speed of light is extremely small at the 
pulsar surface (e.g., Lominadze, Machabeli & Usov 1983), 

c - ?;ph ~ ciojpjoj^^ )^Tp ~ lO-i^c (7) 

in the frame of the pair plasma, particles of both the 
plasma and the beam cannot be in a resonance with t- 
wavcs that propagate in the pulsar plasma, and resonance 
instabilities (for example, the cyclotron instability) do not 
develop in the neutron star vicinities. 

Besides, the energy density of plasma {Sp ~ m^c^npTp) 
is very small in comparison with the energy density of the 
magnetic field (e^ = B'^/Sir), 

ep/e, ~ 2(u;pju;3j% ~ 10-13 (8) 

Therefore, hydrodynamic instabilities are also suppressed 
near the pulsar surface by the strong magnetic field. 

In the near zone of the pulsar magnetosphere (r < 
c/Q), we have Up oc B oc r~^. The ratio i^p/oj^ in- 
creases with increase of the distance from the pulsar, 
u)p/u)^ oc r^/^. Near the light cylinders of pulsars the 
ratio Wp/Wjj is of the order of unity, and the outflowing 
plasma may be unstable with respect to excitation of t- 
waves (e.g., Machabeli & Usov 1979, 1989; Lominadze et 
al. 1983; Lyutikov, Machabeh & Blandford 1999 and ref- 
erences therein). At the distance from the pulsar in the 
range R < r c/Q only the two-stream instability may 
be developed. 

Observations of pulsars suggest that the radio emission 
regions arc far inside the light cylinder (e.g., Gil & Kijak 
1993; Kramer et al. 1994; Kijak & Gil 1997), 

r-radio ^ 50R ~ (0.02 - 0.04) c/0 «: c/Q , (9) 

Therefore, the development of the two-stream instability 

in the magnetospheres of pulsars is the most plausible rea- 
son of their extremely intense radio emission. 



I f(p) 



a. 




Fig. 2. The momentum distribution function of particles 
for a plasma-beam system: (a) prior and (b) after the de- 
velopment of the two-stream instability. 

4. Two-stream instability 

The two-stream instability may occur only if (1) there is a 

"hump" on the momentum distribution function f{p), i.e. 
a region having df/dp > 0, and (2) in the frame where 
plasma is at rest the mean momentum of the "hump" par- 
ticles (pb) is large enough in comparison with the mean 
momentum of the plasma particles (the Penrose criteria). 
Particles, which are in the inside slope of the "hump" 
where df/dp > 0, produce longitudinal (l) waves with the 
frequency co = (k, v) (the Cherenkov resonance), where k 
is the wave vector (see Fig. 1). These particles thus lose 
their energy. As a result the inside slope of the "hump" 
begins to move into the region of even lower velocities, and 
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the distribution function of the " hump" particles becomes 
"plateau-shaped" (see, Fig. 1). The further generation of 
Z-waves is cut off. The plateau formation and the cut-off 
in the plasma wave generation is called the quasi-linear 
relaxation. The energy density of Z-waves generated in the 
process of quasi-linear relaxation is of the order of the 
energy density of the "hump" particles. 

A strong magnetic field has no an influence on the 
development of the two-stream instability in the one- 
dimensional plasma outflowing from the polar cap regions 
of pulsars, and therefore, this instability may occur every- 
where in the pulsar magnetosphere if the mentioned above 
conditions for its development arc realized. Since there are 
several components of particles in the pulsar plasma that 
move with different speeds, the idea that the two-stream 
instability is responsible for the coherent radio emission 
of pulsars is very natural. 

The two-stream instability could be developed due to 
the following three reasons. The first one is the interaction 
between the beam of primary particles and the secondary 
e+e" plasma (Ruderman & Suthermand 1975). The sec- 
ond one is the interaction between electrons and positrons 
of the secondary plasma itself (Cheng & Ruderman 1977). 
At last, the two-stream instability may be triggered in the 
inhomogcncous pulsar plasma when the outflowing plasma 
clouds disperse and overlap each other (Usov 1987; Ursov 
& Usov 1988; Asseo & Melikidze 1998). 




Fig. 3. Charge flow separation of the pulsar plasma mov- 
ing outward along the curved magnetic field lines; /x is the 
magnetic dipole moment. 



4.1. The interaction between the primary beam and the 

secondary plasma 

The ultrarelativistic beam of primary particles has the 
distribution function that obeys the conditions that are 
necessary for the two-stream instability development (sec 
Fig. 1). Therefore, the generation of radio emission of pul- 
sars was initially associated with the ultrarelativistic beam 
of primary particles (Ruderman & Sutherland 1975). How- 
ever, further more precise estimates of the instability 
growth rate have shown that it is impossible. Indeed, the 
characteristic time for development of the two-stream in- 
stability in the pulsar frame is (Benford & Buschauer 1977; 
Egorenkov, Lominadze & Mamradze 1983) 



Ti ~ (np/nfc)i/3r(,ry2^-i . 



(10) 



From equations (l)-(4), (5) and (10), for typical pul- 
sars the value of Tj is about 10~''(r/i?)^/^ s at the distance 
r from the neutron star. For the instability development, 
it is necessary that Ti is, at least, a few times less than 
To = r/c ~ 3 X 10~^ (r/i?) s, the characteristic time of the 
plasma outflow from the distance r. However, the reverse 
relation between Ti and tq holds. Hence, the two-stream in- 
stability that, in principal, could be caused by the primary 
beam because its distribution function forms a "hump" 
with df /dp > 0, in fact, has not enough time to be devel- 
oped in the pulsar magnetospheres. This is because both 



the beam density is too low and the Lorentz-factor of the 
beam is too high. 

4-2. The interaction between the electrons and the 
positrons of the secondary plasma 

The secondary e+e~ plasma produced by photons ini- 
tially has identical distribution functions in its compo- 
nents (electrons and positrons). In this case, the pair 
plasma is neutral, and the charge density pe is determined 
by the primary beam and equals to p^'' = encj- In the 
process of the plasma outflow the charge density pe is 
(xn(x B while p'^^^ is oc B cos x, where x is the angle be- 
tween ft and B. The magnetic field lines are curved, and 
therefore, the angle x varies with the distance from the 
pulsar. This leads to a deviation of pe from and gen- 
eration of the electric field along the magnetic field (see 
Fig. 3). This field accelerates one of the components of 
the secondary plasma and decelerates the other, and the 
distribution functions of these components shift to each 
other. The difference between the mean velocities of the 
electrons and the positrons is (Cheng & Ruderman 1977) 



\v+-V-\ ~ — 



r2 B 

■ B)o 



c, 



(11) 



where (H • B)o is the fi and B product taken at the region 
of pair creation. 




Far the neutron star (r R), equations (3) and (11) yield 

|r+ - r_ 



to r„ 



so that 

Tmax ^ Tjniji ^ 1. Then over the 



c Tb 



where 



(12) 



(13) 



Using a rclativistic generalization of the Penrose crite- 
ria (Buschauer & Benford 1977; Krall & Trivelpiece 1986) 
and equation (12), the condition for development of the 
two-stream instability may be written as 



Tj, > 2 X 102(c/rO)i/3. 



(14) 



Prom equation (14) we can see that the two-stream in- 
stability, which develops because of the relative motion of 
the secondary plasma components, may occur only rather 
near the light cylinder of pulsars. Besides, the magnetic 
field near the pulsar surface must be close to the dipolc 
field when the Lorentz-factor of the secondary particles is 
equal to its maximum (see, equation (1)). 

4-3. The interaction between clouds of the secondary 

plasma 

Here the nonstationarity of plasma ejection from the po- 
lar gap regions is modeled as follows. Plasma particles 

(electrons and positrons) are ejected along the z-axis from 
z = during the time T with Lorentz-factors ranging from 



time T no particles are ejected, after that everything is 
repeated. In other words, we assume that the secondary 
e~^e~ plasma is ejected in the form of separate isolated 
clouds with the length cT along the 2;-axis, and the space 
between these clouds is cr. 

At z = 0, the distribution of ejected particles over 
Lorentz-factors is described by the function /o(r) deter- 
mined from the condition that the number of particles dup 
per volume dV with energies from rUeC^r to mgC^ (F-t-rfF) 
is equal to 



dnp = UpdV fo{T)dr . 



(15) 



The distribution function thus determined is normalized 
to unity, i.e.. 



fo{r)dT 



(16) 



The pulsar plasma can be regarded as collisionless. There- 
fore, the spatial and energy distribution of outflowing par- 
ticles at z > may be found from simple kinematic rela- 
tions. 

To study the dispersion of plasma clouds and their 
mutual overlapping it is sufficient to consider the outflow 
of two successively ejected clouds. It is expedient to draw 
the distribution of particles in the z- and F-coordinates 
(see Fig. 4). In our case the initial time {t = 0) is the 
moment when the ejection of particles of the second cloud 
is just finished. 
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f{r,z,t) 



fo(r) 



Ta r4 



Fig. 5. Overlapping of plasma clouds a,t t > to- The re- 
gions with particles are shaded. In the overlapping region, 
z meets the condition Vmin(T + t + t) < z < fmax(T' + t). 
The distribution function of particles in the overlapping 
region at the distance z from the pulsar surface is shown 
in the bottom. 



At the time 



to 



,(T + t)- 



mm ' 



(17) 



determined from the equation 

Wmin(T + T + to)= W(T + to) , (18) 

the fast particles of the second cloud catch up with the 
slow particles of the preceding cloud, and mutual overlap- 
ping of the clouds begins, where 

^^min=[l-(l/rmm)']'/^C, (19) 
?;max= [l-(l/rmax)^]'/^C. (20) 

The width of the overlapping region a,t t > to is 

(T + t) {T + T + t). (21) 

The distribution function of particles in the overlap- 
ping region is (see Fig. 5) 



f{r,z,t) 



0, r<ri, 

/o(r), ri<r<r2, 

, r2 < Fa , 

/o(F), r3<r<r4, 

U , r4 < F , 



(22) 



where 

Ti = [l-{vi/cf]-^/\ i = 1,2,3,4, 
z 



Vi = v„ 



V3 



i 



T + T + t 

'^4 — '^max* 



(23) 

(24) 

(25) 



In the overlapping region, there are, in fact, two 
streams of slow and fast particles (see Fig. 5). This may 
lead to development of the two-stream instability. A sim- 
ple approximation where plasma is one-dimensional and 
homogeneous may be employed for consideration of this 
instability in the overlapping regions of plasma clouds in 
the magnetospheres of pulsars (Usov & Ursov 1988). In 
this approximation, the dispersion equation for /-waves 
with k II V II B is 



Aire'^np f f{T,z,t)dT 



nif 



(26) 



At the initial stage of clouds overlapping, when the 
width of the overlapping region is much smaller than the 
width of the clouds, along the z-axis, A <C cT, the func- 
tion /(F, z, t) describes two groups of practically monoen- 
ergetic particles at velocities V2 and V4. Taking this into 
account, equation (26) may be written as 



/3 







where 



F3 



r3 



(27) 



(28) 



From equation (27), the maximum growth rate of l- 
waves is (Ursov & Usov 1988) 



Imwmax = -^a^/^ ' ^ ^ 



and occurs for the wave vector 



vl/2. 



Up 



Va - V2 



(29) 



(30) 



that corresponds to a resonance between plasma oscilla- 
tions of the slow (Fi < F < F2) flow and high-energy 
particles of the fast (F3 < F < F4) flow. 

Since for the slow and fast flows the F ranges of plasma 
particles are narrow for slightly overlapping clouds (A <C 
cr), the function /o(F) in equations (28) can be factored 
out the integral sigh and then to integrate: 



a : 



/o(ri 



F2 
^ 2 



V2 — V\ 



/o(ri) 



Ai 

Cto 



/o(rmm) (31) 
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_ /o(r4) rj - t.2 - 1>3 , . . , A2 

P - ^ p2p2 - Jo\>- i) - ^•'0(J- maxj, (J^j 

where (see Fig. 5) 

Ai = - VminCT + r + , A2 = t;max(T + - ^; (33) 
in the overlapping region, 

{T + t). (34) 

Substituting equations (31) and (32) into equation 
(29), we obtain 



V3 



ImWmax = ^^Wp[/o(r„iin)]^^®[/o(r: 



11/3 



A- A2 

Cto 



1/6 /A2^'/' 



(35) 



The growth rate (35) depends on the distanee A2 to 
the leading front of the cloud overlapping and reaches the 
maximum 



(Ima;„,ax)max = ^[/o(r„,i„)]l/'^[/o(r,nax)]'/' 



X (A/cto)V2 



(36) 



at A2 = (2/3)A. 

The distribution function of plasma particles in the 
magnetospheres of pulsars may be roughly approximated 
by the following law (Arons 1981) 



f /p\ ) ^^'^ rjnin < r < r* , 

JoK^ l-S A p-3/2 „ „ „ 

-^i-* max t 



where A is the constant determined from the normaliza- 
tion condition (16) and written as 



^1/2 



A = 



1/2 



-3 



(38) 



Substituting equation (37) into equation (36) and tak- 
ing into account that do — 2Lr^ijj, we have the charac- 
teristic time of the instability development 



Ti — (^imWmaxJmax ~ . . AMI \ ^ 



(39) 



where L = cr is the space between the plasma clouds. 

An instability manages to develop if the lifetime of 
the unstable system exceeds the characteristic time of the 
instability development by more than an order. In our 
case, the unstable system is the overlapping region, and 
its lifetime is A/c. With a similar criterion, lOri = Aj/c, 



wc get the value of A at which the two-stream instability 
is developing to be given by 



(40) 



We now estimate the quantities that characterize the 
development of the two-stream instability in the magne- 
tospheres of pulsars. The following parameters of the sec- 
ondary e'^e~ plasma may be assumed: r„iiji ~ 5, ~ 10^, 
Fmax ^ 10^ (Arons 1981) and L = cr ~ cT ~ i? ~ 10^ 
cm (Usov 1987). Prom equations (3)-(5), (38) and (40), 
we have A^ ~ 0.8 x 10^ cm. Hence, the development of 
the two-stream instability already occurs when the cloud 
overlapping is very small (A, ^ L). Further overlapping 
of the plasma clouds (A > A.;) could be considered only 
beyond the scheme of our linear approximation. 

The distance from the neutron star to the region where 
the two-stream instability occurs is 



n ^ 2crr4i„ = 2Lr4i„ 



(41) 



or numerically ~ 5x 10^ cm for Lx 10^ cm and Fmin — 5. 

The value of is consistent with the distance from the 
neutron stars to the radio emission regions (see equation 
(9)). 

The existence of several scales of the outflowing plasma 
modulation (e.g., Li ~ 0.3it! and L2 — 2Rc in the 
Ruderman-Sutherland model) may lead to the situation 
that in the magnetospheres of pulsars there are several 
regions where the two-stream instability develops and the 
generation of coherent radio emission occurs. Prom the 
observational data on radio emission of pulsars it follows 
that such a situation is possible for some pulsars (e.g., Gil 
1985). 



(37) 5. Discussion 



The mechanism of coherent radio emission of pulsars is 
closely connected with the location of the radio emitting 
regions in the pulsar magnetospheres. If it is confirmed 
that the distances from the neutron stars to the radio 
emitting regions are <C c/il (Gil & Kijak 1993; Rankin 
1993; Kramer et al. 1994; Kijak & Gil 1997), i.e., these 
regions are far inside the light cylinder, then most prob- 
ably, the two-stream instability developed in the outflow- 
ing strongly nonhomogeneous plasma (Usov 1987; Ursov 
& Usov 1988; Asseo & Melikidze 1998) is responsible for 
the generation of the pulsar radio emission. It is worth 
noting that if the process of pair creation near the pulsar 
surface in strongly nonstationary indeed, the development 
of this instability in the magnetospheres of pulsars is al- 
most inevitable. 

Longitudinal (nonescaping) Langmuir waves are gen- 
erated in the process of development of the two-stream 
instability. Then, these waves may be converted into t- 
waves, that can escape from the pulsar magnetosphere, 
by means of the following nonlinear effects: 
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Two-Stream 
instability 


(Ursov& Usov 1988) 


L-waves 


(Asseo, Pelletier & Sol 1990) ^ 


t-waves 





Transformation induced 
by the velocity shear 
(Machabeliet al. 1997) 



L-waves 



Induced scattering 
(Lyubarskii 1996) 



t-waves 



Modulational instability 
and formation of bunches 
(Melikidze et al. 1984; 
2000; Lesch et al. 1994; 
Gil et al. 2000, 2002) 
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